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A METHOD FOR SIMULATING THE ATMOSPHERIC ENTRY OF LONG-RANGE 

BALLISTIC MISSILES 1 


By A. J. Eggebs, Jr. 


SUMMARY 

It is demonstrated with the aid of similitude arguments that a 
model launched from a hypervelocity gun upstream through a 
special supersonic nozzle should experience aerodynamic healing 
and resulting thermal stresses like those encountered by a long- 
range ballistic missile entering the earth’s atmosphere. This 
demonstration hinges on the requirements that model and missile 
be geometrically similar and made of the same material, and 
that they have the same flight speed and Reynolds number 
(based on conditions just outside the boundary layer ) at corre- 
sponding points in their trajectories. The hypervelocity gun 
provides the model with the required initial speed, while the 
nozzle scales the atmosphere, in terms of density variation, to 
provide the model with required speeds and Reynolds numbers 
over its entire trajectory. Since both the motion and aerody- 
namic heating of a missile tend to be simulated in the model 
tests, this combination of hypervelocity gun and supersonic 
nozzle is termed an atmosphere entry simulator. 

INTRODUCTION 

The aerodynamic heating of a long-range ballistic missile 
entering the earth’s atmosphere poses problems of so serious 
a magnitude that the success or failure of the missile may 
well hinge on their solution. The extent to which these 
problems are solved can be fully determined only by flight 
tests, but such tests are, per se, very difficult, time consum- 
ing, and costly. It is appropriate to inquire, therefore, if a 
method can be devised for simulating with relatively simple 
equipment on the ground, the aerodynamic heating and re- 
sulting thermal stresses in a ballistic missile. This question 
forms the starting point of the present paper which under- 
takes first to establish conditions of simulation and then to 
determine a practical method of simulation. In the course 
of events it will be discovered that simulation of heating goes 
hand in hand with simulation of motion of a missile; accord- 
ingly, the subject of this paper is, more generally, an atmos- 
phere entry simulator. 

SYMBOLS 

reference area for drag evaluation, sq ft 
rate constant (see eq. (Al)), sec -1 
rate constant (see eq. (A2)), cu ft/sec 
drag coefficient, dimensionless 
local skin-friction coefficient based on conditions 
just outside of boundary layer, dimensionless 

1 Supersedes NACA Research Memorandum A66116 by A. J. Eggers, Jr., 1S56. 
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Of equivalent skin-friction coefficient, 

~ J O fl dS, dimensionless 

D ' typical model or missile dimension, ft 

E Young’s modulus, lb/sq ft 

H heat transferred per unit area, ft-lb/f t 2 

K diffusivity, sq ft/sec 

k thermal conductivity, ft-lb/sec ft a (°B./ft) 

L length of nozzle, ft 

M flight Mach number, dimensionless 

m mass of missile or model, slugs 

N concentration of particles (i. e., number of parti- 

cles per unit volume), ft -3 
n distance normal to boundary, ft 

h reduced normal distance, n/D, dimensionless 

P 0 nozzle reservoir pressure, lb/sq ft 

Q total heat transferred, ft-lb 

r radius of curvature of body surface at stagnation 

point, ft 

S surface area, sq ft 

T temperature, °R 

T 0 nozzle reservoir temperature, °F 

t time, sec 
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reduced time, jjy sec/sq ft 

velocity, ft/sec 
altitude, ft 

rectangular coordinates, ft (except when appear- 
ing as subscripts on stresses) 
reduced coordinates (jc=x/D, y=y\D, z=z/D), 
dimensionless 

coefficient of thermal expansion, ft/ft- 0 R. 
constant in density-altitude relation, ft -1 (see eq. 

(U) 

typical flow deflection angle, deg 
angle of flight path with respect to horizontal, 
deg 

Poisson’s ratio, dimensionless 
air density, slugs/cu ft 

reference density, slugs/cu ft (reservoir density 
in nozzle) 

tensile or compressive stresses, lb/sq ft (see 
sketch) 

shear stresses, lb/sq ft (see sketch and note 
rxy =z Vy X , Txj—T a , and Tpz— r*j) 
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direction cosines of normal to a boundary, dimen- 
sionless 

SUBSCRIPTS 

ao average conditions 

e conditions at “entrance” to earth’s atmosphere 

or simulator 

l local conditions just outside of boundary layer 

mi missile 

mo model 

s stagnation conditions 

ANALYSIS 

In order to clarify the basic requirements of a simulator, 
this study is initiated with a review of some of the more 
important causes and effects of aerodynamic heatmg of long- 
range ballistic missiles. Following this review, the conditions 
of simulation will be set forth in a mathematical form de- 
signed to facilitate the choice of a practical simulator. 

MOTION AND HEATING OF LONG-BANGE BALLISTIC MISSILES 

It is a fundamental characteristic of ballistic missiles 
(see ref. 1) that speed begets range, with the result that 
hypervelocities in excess of 10,000 feet per second are required 
in order to obtain long ranges in excess of 1,000 miles. A 
long-range ballistic missile first attains hypervelocities near 
the end of powered flight. In this phase of flight and through- 
out the large majority of unpowered flight the vehicle should 
pass more or less unimpeded through the rarefied upper 
atmosphere of the earth, corresponding to altitudes in excess 
of several hundred thousand feet. 1 2 Its trajectory terminates, 
however, with a very rapid descent through all or part of 
the earth’s relatively dense lower atmosphere. In this phase 
of flight, termed the atmospheric entry, retardation and 
severe aerodynamic heating of the missile can almost cer- 
tainly be expected to occur (see ref. 2) . 

Retardation during atmospheric entry is caused by the 
combined action of pressure and viscous forces, while aero- 
dynamic heating stems in the main from work done by 
viscous forces. In both cases it is aerodynamic rather than 
gravity forces which play the predominant role, with the 
result that motion and heating of the missile emerge as 
closely related phenomena. Thus changes in missile shape 
which affect motion will also affect heating. This fact can, 
as discussed in reference 2, be exploited in the design of 
missiles with reduced aerodynamic heating. The potential 
for excessive heating remains, however, as an unavoidable 
property of the long-range ballistic missile which enters the 
earth’s atmosphere at hypervelocity. 

Excessive heating can have several effects. First and 
perhaps most serious of these effects is the development of 
high thermal stresses in the structure of the missile. These 
stresses tend, for example, to far overshadow the pressure 
induced stresses. In addition, there is the natural weakening 
of the missile material at high temperatures, so that struc- 
tural failure may occur during atmospheric entry. There is, 
of course, the further possibility that intense aerodynamic 

1 It Is presumed tliat the rnlsjQe is, unlike the usual meteor, so large (say of the order of feet 

In typical dimension) that free molecule phenomena play a minor role in Its motion and heating. 


heating will, as with meteors, cause burning and the ultimato 
destruction of the missile. ' , 

These, then, are some of the important causes and effects 
of aerodynamic heating of a long-range ballistic missilo 
entering the atmosphere. They suggest that an atmosphere 
entry simulator might logically include means of simulating 
missile velocity, missile configuration (e. g., shape), and the 
lower portion of the earth’s' atmosphere. We will proceed 
from this suggestion to formulate our conditions of simula- 
tion. ; 

CONDITIONS OF SIMULATION 

The analysis to follow tacitly presumes the validity of 
many time-tested assumptions of aerodynamics, thermo- 
dynamics, and solid mechanics. In addition, however, it is 
predicated on the assumptions that during atmospheric 
entry (1) radiation has a secondary effect on missile heating, 
(2) gravity has a secondary effect on missile motion, (3) the 
flow-field freeze principle of Oswatitsch (ref. 3) holds for the 
missile, and (4) the thermal properties (e. g., specific heats) 
and transport properties (e. g., thermal conductivity and 
viscosity) of air are functions of temperature only. The 
first two assumptions are suggested by the calculations of 
references 1 and 2, 3 while the third assumption hinges essen- 
tially on the requirement that the square of the hypersonic 
similarity parameter for a missile be large compared to 1 
during entry (i. e., M 2 sin 2 5»1). For missiles of normal 
slenderness, the hyperspeed of entry tends to insure the 
satisfying of this requirement with the result that flight 
Mach number loses its significance as an important similarity 
parameter.* The last assumption is restrictive only in the < 
event air tends to dissociate (or possibly ionize) and it will 
therefore be treated in a discussion of those phenomena 
later in the paper (see section on “Performance Limitations”) . 

Now it is convenient in discussing similitude to imagine a 
model counterpart to the missile and a test chamber counter- 
part to the atmosphere. Furthermore it is permissible for 
our purposes to proceed from the simplified equations of 
reference 2 for the convective heating of ballistic missiles. 
Thus for an isothermal atmosphere (which closely approxi- 
mates the earth’s lower atmosphere, see ref. 2) 



and there follows: 

(a) Heat absorbed per unit mass at altitude y' , 

(2> 

where 

i £pPo£_ Ml 

yi = ytg pm sin 6, \6) 


3 The first assumption Is, of course, best suited to ‘'relatively light" missiles whloh aro de- 
signed on the so-called “heat sink" principle, or, more generally, which are deslgnod to main- 
tain rel ativoly cool surfaces. 

* It Is demonstrated in reference 3 that Mach number and stream angle In tho disturbed 
flow are independent of flight Mach number, provided AP sin* J»l, and provided the air 
behaves Ideally. In the event nonideal behavior, like changes In the specific heats, occurs 
as a result of high temperatures In the disturbed flow, then flight velocity tends to replace 
flight Mach number as the Important Index of motion since it is through this velocity (via 
kinetlo energy) that high disturbed air temperatures are created. 
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(b) Average rate of heat transfer per unit area 

3 Ojy p<yA 


JT7 P ' n 17 3 -- DP — 

- W Po y 2 sin 0, 




dt 4 

(c) Rate of heat transfer to stagnation point 

itt I /Sy' 3 , 

i~i==const F/e 2 e~^ m sln *• 


(4) 


(5) 


According to equation (1) the density of the air in the test 
chamber must vary exponentially with the distance corre- 
sponding to altitude in the atmosphere. The more general 
implication is, of course, that the test chamber must duplicate 
variations of p/p 0 in the atmosphere, whatever they may be, 
although the absolute magnitudes of p and p 0 may be quite 
different from those in the atmosphere. The static tempera- 
ture of the air in the test chamber is, as in the case of the 
atmosphere, presumed to be small- by comparison to missile 
recovery temperature (see ref. 2 in connection with this point 
as it relates to the derivation of eqs. (2) through (5)). 

It will be stipulated now that model and missile be geome- 
trically similar in structure and configuration, and made of 
the same material. Furthermore, the condition is imposed 
that the model enter the test chamber at the same speed and 
temperature as the missile enters the atmosphere. Finally, 
it is required that model and missile have the same Reynolds 
numbers (based on local conditions outside the boundary 
layer) at corresponding points fiy' in their trajectories. By 
corresponding points it is meant where the product fiy' is the 
same for model and missile. It should be recognized, of 
course, that, in general, fi and y' will individually be grossly 
different for model and missile. 

It follows from these requirements and equation (2) that 
the heat transfer per unit mass Q/m to model and missile 
will be the same at corresponding points fiy' provided V is 
the same, since 0/ and O/S/CoA are the same. But from 
equation (3) the velocity V will be the same at corresponding 
fiy' only if C D p 0 A/fim sin 6, is the same. If the subscript mo 
refers to model and mi to missile, then the last provision may 
bo written. 

(y'/am d,) ma , R . 

(p a D)mc \D m J (y'/sm d t ) ni 


But model and missile Reynolds numbers, velocities, and 
disturbed air temperatures 6 are the same, hence 

( pJD)mt={pJD) m o (7) 

and equation (6) may be written 



* This observation with regard to disturbed air temperatures la easily verified by considering 
flow near the surface of a missile with a stagnation point. Assuming for simplicity that air 
In the disturbed flow behaves Ideally, we have at the stagnation point T.ra VyiC, since 
M*» 1. Then the temperature of the air Just outside the boundary layer Is given by the 

y l / v—i \-i 

expression Tr— ^,-1 1 | 2 M fj • Hence If Vand Mi are the same for model and missile. 

then Tx and Vi are the same, Independent of ambient air temperature. We are assured of 
euunl Mi a by the freme principle (see footnote 4). 


This expression fixes the length L of the test chamber in 
terms of the portion of the atmosphere to be simulated 
therein and the ratio of model to missile size.® If equation 
(8) and the previously set forth requirements are satisfied, 
then model and missile should experience equal heat, transfer 
per unit mass, and hence equal average temperature rise at 
corresponding points in their trajectories. These quantities 
are significant, of course, because they tend to determine 
whether a missile will melt or perhaps bum during flight. 

The next question is how do the heat-transfer rates com- 
pare in the case of model and missile? It is easily deduced 
from equations (4) and (5) and the conditions for equal 
heat transfer per unit mass that 


(dE a fi 

\ D mi , 

fdH m \ 

(9) 

V dt , 

Lo D„o 

l dt ) mt 

(dE£ 

\ D m i 1 

fdH,\ 

(10) 

V At y 

) ~d ' 

■mo -Lsmo 

V dt ) ni 


at corresponding points fiy’. That is to say, the average 
and stagnation point heat-transfer rates are higher for the 
model in proportion to the ratio of missile to model size. 
But perhaps the foremost importance of heat-transfer rates 
is, as discussed earlier, in how they influence thermal stresses 
in the missile structure and, for example, lead to ablation 
of surface material. Evidently, then, it would be most 
desirable if equations (9) and (10) implied equal thermal 
stresses in model and missile. This possibility is easily 
checked using modified equilibrium thermal-stress equations 
for an unrestrained isotropic elastic body (see ref. (4) and 
sketch). 


z 



8 Since model and missile velocities are the same at corresponding points in their trajectories. 
It follows from equation (8) that the time of flight in the test chamber is reduced below that of 
atmospheric entry by the ratio 
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dfrr . dr CT . dr„ aE bT =n ' 
dx by bz 1— 2v dx 

Z><r v ■ dr CT dry, aE dT n n 

dy^ dx "*■ dz 1 — 2v by | 

ba~i . br rt br„ aE bT „ 
dz dx ”* by 1— 2v bz 

where the stresses are not true thermal stresses, but rather 
they are the stresses produced by the “body forces,” 
~aEbT , 

1 — 2r bx’ etC ‘ 

The corresponding compatibility conditions are 

7V *+tb (jh v ‘ T + 2 w) 

- (12a) 

1 bV aE &T 

T “ l+p bybz 1— 2v bybz 

w , i av -_o ^ yy 

djds ^ ! — 2 p bxbi 


where 


ui , i by «g b*r 

Tar ^l-f r bxby 1— 2v bxby 


v dx’^dp^dz* 

^Oi+O-J+O-, 


and the boundary conditions are 

>i , aET ' 

&x + i-^xyV ”1 Trj i — ^ 9^ s 

. . f. aET 

o»»J+t , m ?+t j ,?=j = ^ ij - 

. v , g-ST 

°~»r+ + ?vi0 = j- — g - r 

The"lieat-flow equation is (see ref. (5)) 

b*!T b a T b 3 T 1 bT 
V bx s 1 bp + bz s K bt 


with boundary conditions 

(T aa u) in«!oj=known 


, bT r.dH 
k =D-j— 

071 0 uUr Clt 

surface 

=o 

071 inner 
surface 


and (probably) 


Now the thermal stresses are given by superposing the 
"hydrostatic” pressure aETJl—2v on solutions to equations 
(11) consistent with the compatibility and boundary con- 
ditions (eqs. (12)) and the solution to the heat-flow equation 
(13) along with its boundary conditions (eqs. (14)). But if 
equations (9) and (10) along with the requirements necessary 
for their development are satisfied, then equations (11) 
through (14) (and the “hydrostatic” pressure) are mathe- 
matically identical for model and missile. Hence, the ther- 
mal stresses must be the same in model and missile at corre- 
sponding f3y' in their trajectories. 

Necessary conditions for simulating aerodynamic heating 
and resulting thermal stresses in a ballistic missile are then, 
according to this analysis, as follows. Fust, model and 
missile must be geometrically similar and made of the same 
material. In addition, they must have the same flight 
speeds and Reynolds numbers (based on local conditions 
outside the boundary layer) at corresponding points in their 
trajectories. Finally, in order to meet these conditions and 
to insure equal heating of model and missile, the test chamber 
must contain air at relatively low temperature and with 
variations in p/p 0 along the model flight path equal to those 
in the atmosphere along the missile flight path. 

With this knowledge we are in a position to consider the 
practical problems of simulating atmospheric entry. 

" A PRACTICAL ATMOSPHERE ENTRY SIMULATOR 

It is appropriate to determine first how to provide a model 
with the correct initial hypervelocity required for simula- 
tion. For this purpose it is suggested that a hypervelocity 
gun can be employed. Several such guns have been developed 
(see e. g., ref. 6) to launch small models at speeds above 
15,000 feet per second. These speeds correspond to those of 
ballistic missiles (see ref. 1) with ranges of the order of 2,000 
miles and greater. Both speed and range fall, then, into 
the categories of interest in this paper. 

The next question is how to provide a model test chamber 
which simulates the lower portion of the earth’s atmosphere. 
For this purpose it is suggested that a special supersonic 
nozzle can be used to advantage. To illustrate, it was found 
in reference 2 (see fig. 4 therein) that the major part of the 
aerodynamic heating of a ballistic missile entering the atmos- 
phere occurs over a range of about 100,000 feet in altitude. 
In this altitude range p/p 0 varies by a factor of about 10~ s . 
A corresponding variation in density can be obtained be- 
tween the settling chamber and exit section of a Mach 
number 5 supersonic nozzle (see ref. 7). Imagine, then, a 
hypervelocity gun positioned to launch a model upstream 
along the axis of such a nozzle. The nozzle is designed to 
provide an essentially exponential variation in density along 
its axis, the density decreasing from a very high value in 
the settling chamber to a very low value at the exit. Accord- 
ingly, a model proceeding upstream through the nozzle en- 
counters an increasingly dense atmosphere like that pre- 
sented by the earth to a descending ballistic missile. Now, 
to be sure, unlike atmospheric air, the air in the nozzle is in 
motion. However, the air velocity is small by comparison 
to the hypervelocity of the model, and therefore this differ- 
ence between nozzle air and atmospheric air should mar 
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only slightly the function of the nozzle as a test chamber. 7 

The combination of hypervelocity gun and supersonic 
nozzle merits attention, then, as an atmosphere entry simu- 
lator. An example simulator of this type will therefore be 
considered next. 

EXAMPLE SIMULATOE 

Consistent with the previous discussion, a Mach number 5 
nozzle is chosen to simulate the earth’s atmosphere over 100,- 
000 feet of altitude. Now let us assume that ballistic missiles 
with up to 4,000 miles range are to be studied with this simu- 
lator. The corresponding range of atmospheric entrance an- 
gles 6 e is from 45° down to about 30° (see ref. 1). The length 
of the nozzle is fixed according to equation (8) by the maxi- 
mum values of (y//sin 6,) mi and D mo /D ml . This value of 

(y//sin 0 e ) mt is, from the above specifications, ' „ = 

200.000 feet. The maximum value of D mc ID mi will be dic- 
tated by the size of the largest model which can be launched 
by the special gun available and by the size of the smallest 
missile to be simulated. For the purpose of this discussion 
it suffices to observe that the model size will probably be of 
the order of a fraction of an inch, while the missile size will 
probably be of the order of several feet. It follows that the 
maximum value of D m0 /D m{ should be of the order of 10 -2 . 
In this event we have from equation (8) that the length of 
the nozzle is of the order of 20 feet, and from equation (7) 
the nozzle stagnation density is of the order of 100 times sea 
level ah’ density. A missile 3 feet in diameter and weighing 

5.000 pounds would thus be simulated by a model 0.36 inch 
in diameter and weighing 0.005 pound. 

On the basis of these considerations our example simulator 
might appear something like the one shown in figure 1. The 
required nozzle stagnation densities are obtained for running 
times of the order of a second by using a settling chamber 
charged with high pressure air. A wide range of settling- 
chamber pressures and, hence, densities makes it possible to 
vary the range of altitudes simulated. The storage heater 
in the settling chamber maintains the air temperatures above 
those for which liquefaction can occur in the nozzle. The 
nozzle contour will tend to have very small slopes in the 

> It la Interesting to note that tha supersonic noule when used in the proposed manner has 
compensating features; namely, as model velocity decreases, the air velocity decreases and, 
as model recovery temperature decreases, the air static temperature rises up toward atmos- 
pborio air temperature. These features are, of course, favorable to simulating atmospheric 
entry In that they tend to preserve the required similarity between atmospheric air and noirte 
air. In this manner, too, the possibility of correcting for small differences between these 
media is enhanced. 


streamwise direction. Hence it is anticipated that flow in 
the nozzle will not depart radically from the one-dimensional 
type. Air from the nozzle passes into a vacuum tank which 
is of sufficient size to maintain the nozzle compression ratios 
required for supersonic flow during the course of a test. Ob- 
servation windows are located at short intervals along the 
nozzle side walls to permit photographing models in flight 
and determining effects of aerodynamic heating. In this re- 
gard it would quite likely prove desirable to employ a spec- 
trometer to identify the sources of radiation energy emitted 
in the vicinity of models. The simulator would be instru- 
mented for measuring settling-chamber pressure and temper- 
ature and, of course, the time-distance history of models. 
Models would be launched from a hypervelocity gun located 
at or near the end of the nozzle, and they would be arrested 
in a catcher located in the settling chamber. The damage to 
models in the catcher should be negligible for low impacting 
velocities, thereby permitting the study of recovered models 
to determine changes which occurred during atmospheric 
entry. 

We have, then, some idea of how the proposed simulator 
might be employed in practice. It is important, however, 
to be aware of the points of difficulty which may limi t the 
performance of the device as a simulator. 

performance limitations 

First, it should be noted that pressure-induced stresses in 
the model will be higher than those in the missile by the 
ratio DmilDmo. This point may prove troublesome, al- 
though not unduly so if these stresses are very low (as they 
tend to be) in the missile. It should be remembered too 
that by proper model design (e. g., the use of internal 
pressurization or other bracing) this problem can be mini- 
mized. Careful attention must, of course, be given to both 
model and sabot design from the standpoint of minimizing 
stresses in the model during launching. 

Another point of difficulty may be encountered if material 
properties (e. g., yield point) are a significant function of 
time under rapid heating conditions. There is some indica- 
tion (see ref. 8 and papers cited therein) that at the very 
high heating rates of long-range ballistic missiles (corre- 
sponding to temperature-rise rates of the order of 100° F/sec 
and more), the importance of time is small. It is indicated 
too that this remains true, and, more important, that 
material properties remain essentially the same at the much 
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higher heating rates produced in the simulator. 8 Far more 
information is needed, however, before the significance of 
time in the sense of this discussion can be fully assessed. 

A further question which should be raised is whether or not 
a body can distort fast enough with changes in temperature 
to remain in stress-strain equilibrium (i. e., acceleration 
terms negligible) as was assumed in deriving equations (11) 
and (12). This situation would be more serious in the case 
of the model than in the case of the missile. A crude check 
on the matter is easily obtained for a body which behaves 
elastically. To illustrate, the time required for the model 
to make small adjustments in stress-strain equilibrium 
should be of the same order as the time required for an 
elastic wave to travel the length of the model. For atypical 
steel model (D=% in.) in the simulator, this latter time 
would be about 10“° seconds. If we assume that the model 
experiences a total temperature rise of 1500° F while traveling 
through the simulator, then 10 -8 seconds is also the time 
required for only about 1° temperature change in the model 
material. Accordingly, only small adjustments in stress- 
strain equilibrium are evidently required in this time, and 
hence equilibrium should tend to be realized in model as well 
as missile. If time-dependent plastic deformation should 
become significant, then this statement obviously no longer 
holds. 9 Thus, for example, the simulator may not (in view 
of its foreshortened time scale) duplicate more than quali- 
tatively a fracture process (see ref. 9) although it should tend 
to duplicate thermal deformation up to and including the 
beginning of fracture. If fracture occurs on a molecular 
scale, like ablation by sublimation, then the phenomenon 
will tend to be properly reproduced in the simulator provided 
the sublimation rate is proportional to heat- transfer rate. 

There is, too, the possibility of a missile being aerody- 
namically heated to temperatures where it will bum during 
descent through the earth’s atmosphere. The simulator 
should tend to duplicate conditions leading up to this 
phenomenon; however, there is a question as to how well 
burning would be duplicated. The complexity and lack of 
complete understanding of metal burning (see, e. g., ref. 10) 
preclude the possibility of obtaining a quantitative answer 
to this question at the present time. From the qualitative 
viewpoint it is reasonable to expect that the effects of 

• This discussion should not bo construed to mean that material properties are the same 
under conditions of high rates of heating as they are under steady state conditions (see ref. 8) . 

• Thus with a slow process, like creep, model deformation would not simulate missile 
deformation (see, e. g. f ref. 9), but rather. It would be less by the ratio of (.DmiJDmd* due to the 
foreshortened time scale In the simulator. Even In the case of the missile, however, the 
entry time may be so short (of the order of seconds) that creep plays a minor role In deforming 
the vehicle. 


increased partial pressure of oxygen (acting to increase 
burning rate) and the reduced model size (acting to decrease 
the amount of material to be burned) should combine in the 
simulator to compensate for its foreshortened time scale, 
thereby more nearly providing simulation of missile burning. 
Burning should be accurately simulated for cases where 
burning rate is directly proportional to the partial pressure 
of oxygen. 

As a final point, it is appropriate to consider the matter of 
dissociation and association of the oxygen and nitrogen in 
air. The simulator produces essentially the same disturbed 
air temperatures as the missile, entering the atmosphere. 
Accordingly, the potential for dissociation and association is 
duplicated by the simulator. At the present time, however, 
it is felt that these two phenomena obey different rate laws 
(see Appendix). In this event, neither the phenomena nor 
their effects on heating of a ballistic missile can be duplicated 
except by 1 to 1 simulation (i. e., the equivalent of flight tests 
with the full-scale missile). An indication of the possible 
error that this situation may introduce in tests with the pro- 
posed simulator can he obtained from the calculations of 
reference ll' which suggest that the net effect of equilibrium 
dissociation and association in free flight may be to increase 
the rate of heat transfer to a stagnation point by only about 
50 percent. 10 This increase should be well within the design 
safety factor of a missile. 

CONCLUDING REMARKS 

It has been found that an atmosphere entry simulator con- 
sisting of a hypervelocity gun combined with a special super- 
sonic nozzle may be used to study the conditions of extreme 
heat transfer and thermal stress which introduce such serious 
problems in the flight of a long-range ballistic missile. The 
effects of aerodynamic heating on the model can be observed 
with relative ease, and further, the tests can be conducted at 
a cost which is negligible by comparison to that of flight 
tests. Indeed, in the simplest test, the simulator could pro- 
vide with one photograph of a model rather substantial evi- 
dence as to whether or not the corresponding missile would 
remain essentially intact while traversing the atmosphere. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 15, 1955 


n Dissociation and association may nave especially marked effects on stagnation point heat 
transler since these phenomena can strongly affect both the invisdd and viscous flows In this 
region. 
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APPENDIX A 

DISSOCIATION AND ASSOCIATION RATES 


The dissociation of 0% and N 2 is thought (see, e. g., refs. 12 
and 13) to obey the linear law 

^=a(T)N (Al) 

where N is the concentration (or density in particles per unit 
volume) and Oi(T) is the rate constant which usually de- 
pends only upon temperature. On the other hand, the cor- 
responding association process (e. g., O+O-^Oi) is thought 
to obey the second-order law 1 

( —=C 2 (T)N 1 (A2) 

Now disturbed ah 1 velocities are the same in the case of model 
and missile, while disturbed air densities are higher in the 
case of the model by the ratio of D mi jD n0 . It follows then 
that the simulator will tend to duplicate a rate process in 
which the percentage rate of change of concentration of a 
given type of particle is proportional to the concentration of 
that particle, namely, 

j}f-C(T) N (A3) 

But the simulator tends to duplicate temperature and type 
of particle; hence, it should duplicate the “association” rate 
process given by equation (A2) since C(T)=Ot{T). It will 
apparently not, however, duplicate the “dissociation” rate 
process given by equation (Al). 


l Actually there Is some question as to the correctness of either equation (Al) or (A2) (see 
o. g„ ref. (14)) for pure media and there Is the further complication of Impurities which could 
lead, for example, to a third-order rate process. 
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